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Abstract. The field theoretic renormalization group and the operator product 
expansion are applied to the stochastic model of passively advected vector field with 
the most general form of the nonlinear term allowed by the Galilean symmetry. The 
advecting turbulent velocity field is governed by the stochastic Navier-Stokes equation. 
It is shown that the correlation functions of the passive vector field in the inertial 
range exhibit anomalous scaling behaviour. The corresponding anomalous exponents 
are determined by the critical dimensions of tensor composite fields (operators) built 
solely of the passive vector field. They are calculated (including the anisotropic sectors) 
in the leading order of the expansion in y, the exponent entering the correlator 
of the stirring force in the Navier-Stokes equation (one-loop approximation of the 
renormalization group). The anomalous exponents exhibit an hierarchy related to 
the degree of anisotropy: the less is the rank of the tensor operator, the less is its 
dimension. Thus the leading terms, determined by scalar operators, are the same as in 
the isotropic case, in agreement with the Kolmogorov's hypothesis of the local isotropy 
restoration. 



1. Introduction 

In the past two decades, much attention has been attracted by turbulent advection of 
passive scalar fields; see the review paper and references therein. Being of practical 
importance in itself, the problem of passive advection can be viewed as a starting point 
for studying intermittency and anomalous scaling in the fiuid turbulence on the whole 
[2]. Most progress was achieved for the so-called Kraichnan's rapid-change model: for 
the first time, the anomalous exponents were derived on the basis of a microscopic 
model and within controlled approximations [3J. In Kraichnan's model the advecting 
velocity field Vi{x) with x = {t,x.} is modelled by a Gaussian statistics with vanishing 
correlation time and prescribed correlation function (vv) oc 6{t — t')k~'^^^, where k is 
the wave number, d is the dimension of space and ^ is an arbitrary exponent with the 
most realistic (Kolmogorov) value = 4/3. 

The "zero-mode approach," developed in [3], is based on the fact that, due to 
the temporal decorrelation of the Kraichnan's velocity field, some closed differential 
equations can be derived for the equal-time correlation functions of the passive fields. 
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In this sense, the model is equivalent to a certain quantum-mechanical problem and 
appears "exactly solvable." Although such equations cannot be explicitly solved, 
the anomalous exponents can be extracted from the analysis of their infrared (IR) 
asymptotic behaviour; see [1] for a detailed review and the references. 

In [1] and subsequent papers, the field theoretic renormalization group (RG) and 
the operator product expansion (OPE) were applied to Kraichnan's model; see [5] for 
a review and the references. In that approach, the anomalous scaling emerges as a 
consequence of the existence in the corresponding OPE of certain composite fields 
("operators" in the quantum-field terminology) with negative dimensions, which are 
identified with the anomalous exponents. This allows one to construct a systematic 
perturbation expansion for the anomalous exponents and to calculate them up to the 
order [4] and [6]. Besides the calculational efficiency, an important advantage of 
the RG-I-OPE approach is its relative universality: it can also be applied to the case of 
finite correlation time or non-Gaussian advecting fields. 

For passively advected vector fields, any calculation of the exponents for higher- 
order correlations calls for the RG techniques already in the O(^) approximation. 

Owing to the presence of a new stretching term in the dynamic equation, the 
behavior of the passive vector field appears much richer than that of the scalar field: 
"...there is considerably more life in the large-scale transport of vector quantities" (p. 
232 of Ref. [2]). Indeed, passive vector fields reveal anomalous scaling already on the 
level of the pair correlation function O HI |9] . They also exhibit interesting large-scale 
instabilities that can be interpreted as manifestation of the dynamo effect [71 HDl [ll] ■ 
Some special case {A = 0, see below) reveals a close formal resemblance to the NS 
turbulence [I2l[l3]. From the physics viewpoints, passive vector fields can have different 
physical meaning: magnetic field, weak perturbation of the prescribed background fiow, 
concentration or density of the impurity particles with an internal structure. 

In this paper, we study anomalous scaling of a passive vector quantity, advected 
by a non-Gaussian velocity field with finite correlation time, governed by the stirred 
Navier-Stokes (NS) equation. 

The plan of the paper is the following. In sec. [2] we give detailed description 
of the stochastic problem and explain the motivation of our study and its relation 
to previous work. In sec. |3] we give the field theoretic formulation of the model. In 
sec. m we establish its renormalizability, derive the corresponding RG equations and 
present the explicit one-loop results for the renormalization constants and RG functions. 
Possible IR attractive fixed points are discussed in sec. |5l In sec. [6] the operator product 
expansion is employed to establish the anomalous scaling of the correlation functions 
in the inertial-range. The corresponding anomalous exponents are determined by the 
critical dimensions of tensor composite operators built solely of the passive field. The 
practical calculation is performed in the leading (one-loop) approximation; the results 
are presented in sec. [3 Section [8] is reserved for a brief conclusion; in particular, we 
mention an hierarchy demonstrated by the anisotropic contributions. 
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2. Description of the model 

We confine ourselves with the case of transverse (divergence-free) passive 6i{x) and 
advecting vector fields. Then the general advection-diffusion equation has the 

form 

Vte^-Ao{ek^k)vi + ^iV = ^^o^^ei + r]i, Vt = dt + {vkdk), (2.1) 

where Vt is the Lagrangian (Galilean covariant) derivative, V{x) is the pressure, kq is 
the diffusivity, is the Laplace operator and rii{x) is a transverse Gaussian stirring 
force with zero mean and covariance 

{m{x)Vk{x')) = S{t - t') C,u{r/L). (2.2) 

The parameter L is an integral scale related to the stirring, and Cik is a dimensionless 
function with the condition diCik = 0, finite at r = and rapidly decaying for r — )■ oo; 
its precise form is unessential. Due to the transversality conditions diO^ = 0, diVi = 0, 
the pressure can be expressed as the solution of the Poisson equation: 

d^V = (A - 1) diVkdkOi. (2.3) 

Thus the pressure term makes the dynamics (12.11) consistent with the transversality. 
The amplitude factor Aq in front of the "stretching term" {6kdk)vi is not fixed by the 
Galilean symmetry and thus can be arbitrary. Such general model" was introduced 
and studied in refs. [14]-[16]; it can be naturally justified within the so-called multiscale 
techniques. 

From the physics viewpoints most interesting is the special case Aq = 1, where the 
pressure term disappears: it corresponds to magnetohydrodynamic (MHD) turbulence 
^17j . It was studied earlier in numerous papers; see e.g. refs. [3 EJ [9l [HI [19] and 
references therein. 

In earlier studies, the velocity field in (12. ip was usually described by the Kraichnan's 
rapid-change model: Gaussian statistics with vanishing correlation time and prescribed 
power-like correlation function. In this paper, we employ the stochastic NS equation: 

VtVi = uod\-dip + fi, (2.4) 

where is the same Lagrangian derivative, p and fi are the pressure and the transverse 
random force per unit mass. We assume for / a Gaussian distribution with zero mean 
and correlation function 

{Mx)fj{x')) = ^-^^f^ [ rfki^,(k)rf^(A;)exp[ik(x-x')], (2.5) 

['^'^) Jk>m 

where -Pij(k) = Sij — kikj/k"^ is the transverse projector, df{k) is some function of A; = |k| 
and model parameters. The momentum m = 1/L, the reciprocal of the integral scale L 
related to the velocity, provides IR regularization. For simplicity, we do not distinguish 
it from the integral scale related to the scalar noise in (12. 2p . 

The standard RG formalism is applicable to the problem (12. 4p . (12. 5p if the 
correlation function of the random force is chosen in the power form [20] 

df{k) = Dok^-'^-\ (2.6) 
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where > is the positive amphtude factor and the exponent < ?/ < 4 plays the 
role of the RG expansion parameter. The most realistic value of the exponent is y = 4: 
with an appropriate choice of the amplitude, the function (12. 6p for y ^ 4 turns to the 
delta function, df{k) oc S(k), which corresponds to the injection of energy to the system 
owing to interaction with the largest turbulent eddies; for a more detailed justification 
see e.g. [2T|l22]. 



3. Field theoretic formulation 

According to the general theorem (see e.g. [21]), the full-scale stochastic problem (12. ip - 
(12. 6p is equivalent to the field theoretic model of the doubled set of fields $ = {v, v', 6, 6'} 
with the action functional 

s{^) = 5,(v', v) + e'Dge'/2 + e' {-Vt - AoiOkd^H + koO^} e, (3.1) 

where Dq is the correlation function (12. 2p of the random force rji in (12. ip and is the 
action for the problem (I2.4p - (l2.6p : 

5,(v', v) = v'DW/2 + v' {-Vt + i^od^} v, (3.2) 

where is the correlation function (12. 5p of the random force /j. All the integrations 
over X = {t, x} and summations over the vector indices are understood. The auxiliary 
vector fields v', 6' are also transverse, diV^ = diO'^ = 0, which allows to omit the pressure 
terms on the right-hand sides of expressions (13. ip . (13. 2p . as becomes evident after the 
integration by parts. For example. 



dt J dx v'^dip = — J dt J dx p{div[) = 0. 



Of course, this does not mean that the pressure contributions are unimportant: the 
fields v', 9' act as transverse projectors and select the transverse parts of the expressions 
to which they are contracted. 

The part of the coupling constants is played by the three parameters go = Dq/uq, 
Aq and Mo = Ko/t'O) the analog of the inverse Prandtl number in the scalar case. By 
dimension, 

go oc A^, ^0 and Uq oc A°, (3.3) 

where A is the characteristic ultraviolet (UV) momentum scale. Thus the model (13. ip . 
(13. 2p becomes logarithmic (all the coupling constants become dimensionless) at y = 0, 
and the UV divergences manifest themselves as poles in y. 



4. Renormalization and RG equations 

The renormalization and RG analysis of the model (13. ip . (13.20 are similar to that of the 
scalar advection by the NS velocity field [23], and here we discuss them only briefly. 
Dimensional analysis shows that superficial UV divergences can be present only in 
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the 1-irreducible Green functions (v'v), (v'vv), {9' 9) and {O'vO). The corresponding 
counterterms reduce to the forms v'dfV, v'd'^v, v'{vd)v, O'dtO, 6'd'^6, 6'{6d)v and 9'{vd)9. 

The spatial derivative d at the vertices v'{vd)v, 0'{6d)v and 6'{vd)6 in (13. ip can 
be moved, using the integration by parts, onto the auxihary fields v' and 6' . Thus 
any counterterm must include one spatial derivative per each auxiliary field. This 
excludes the counterterms v'dtv and 9'dt9 without a spatial derivative. Then the Galilean 
symmetry excludes the structures v'{vd)v and 6'{vd)6 because they must enter the 
counterterms only in the form of Galilean invariant combinations v'VtV and O'Vtd- 

The remaining three counterterms v'd'^v, O'd^O and 6'{6d)v can be reproduced by 
multiplicative renormalization of the parameters 

z/Q = i^Z^, Ko = kZ^, Ao = AZ_A, go = gfi^Zg, Zg = Z~^; (4.1) 

no renormalization of the fields $ = {v, v', 6, 6'} and the IR scale m is needed. Here u, 
g, K and A are renormalized analogs of the bare parameters vq, go, no and Ao, while the 
reference scale fi is an additional parameter of the renormalized theory. The last relation 
in (14. ip follows from the absence of renormalization of the amplitude Do = goi^o ~ 91^^^^ 
in the first term of the action S^r- The renormalization constants Zi = Zi{g, u. A, d, y) 
absorb all the UV divergences, so that the Green functions are UV finite (that is, finite 
aX y = 0) when expressed in renormalized parameters. 
The corresponding renormalized action has the form 

Sr{<I>) = S.niV, v) + e'Dee'/2 + 6' {-V* - AZ^{eA)v, + t^Z^d'} 9, 
S,r{w',^^) = v'D,v'I2 + v' {-Vt + uZ,d^]v, (4.2) 

where is expressed in renormalized parameters using (14. ip . It differs from the original 
(unrenormalized) action (13. ip . (13. 2p only by the choice of parameters, iSR($,e, /x) = 
iSo($,eo), where Cq is the full set of bare parameters and e are their renormalized 
counterparts. Thus the original G = ($... $) and the renormalized Gr Green functions 
are also related as G{eo, ■ ■ ■) = Gr^c, /i, . . .); the ellipsis stands for the other arguments 
(times/frequencies and coordinates/momenta). We use to denote the differential 
operation fid^ at fixed bare parameters Cq and operate on both sides of the last relation 
with it. This gives the basic RG differential equation: 

{V^ - 7,1), + Pgdg + Mu + I^AdA] = 0. (4.3) 

Here u = n/v and Vg = sdg for any variable s. The RG functions (the /3 functions and 
the anomalous dimensions 7) are defined as 

-fF = V^lnZF (4.4) 

for any quantity F and 

/3g = ^ti9 = ^[-l/ + 37^], 
X = V^u = m[7, - 7«], 

= = -A^A (4.5) 
for completely dimensionless variables (coupling constants). Here is the operation 
at fixed bare parameters and the second relations in (14. 5 p follow from the definitions 
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and the relations ( I4.ip . It remains to note that the differential operator in fl4.3p is 
nothing other than expressed in renormalized variables. 
The one-loop calculation gives: 

u[u + 1)^ 2d{d + 2) y 



where 



Q = {u + 1){3A^ + Ad-2A + d'^ -3)- 2A{A - 1), (4.7) 

Sd = Sd/{27rY and 5*^ = 2'n'^/'^ /V{d/2) is the surface area of the unit sphere in d- 
dimensional space. Of course, due to the passivity of the field 6, the constant Zy is 
the same as in the model (13.21) : it does not depend on the parameters u and A related 
to the passive field. It is noteworthy that the expression (14. 7p for Q simplifies for the 
aforementioned special values of A: 

Q = {u + l){d - l){d + 2) for ^ = 1, 

g = (M + l)(rf^-3) for^ = 0. (4.8) 

From (14. 6 p we obtain the following explicit one-loop expressions for the anomalous 
dimensions: 

7. = ^^,^^^ + 0(^?2), 7^ = 0(/), 

^^ = ^^^2^^iT-v^+0^3\ (4.9) 

u[u + 1)^ 2d[d + 2) 

In the rapid-change version of our model [H], Z_4 = 1 and 7^1 = identically 
because all nontrivial Feynman diagrams of the 1-irreducible Green function (O'vO) 
contain closed circuits of retarded propagators and therefore vanish. In the present 
case, the absence of the 0{g) term in Z^ and 7^ is a result of the cancellation of 
the (nontrivial!) contributions from the three one- loop diagrams in the 1-irreducible 
Green function {9'v9). For the counterterm 6'{vd)6 such a cancellation is guaranteed 
by the Galilean symmetry (to all orders in g; see the discussion above). For 9'{9d)v 
the cancellation looks accidental and can be explained by a rather simple form of the 
one-loop diagrams: the structures corresponding to the counterterm 9'{vd)9 cancel each 
other due to the Galilean symmetry, while the structures corresponding to 9'{9d)v enter 
all the one-loop diagrams with the same coefficients and cancel out into the bargain. 
This mechanism is not expected to work beyound the one-loop approximation; thus 
nontrivial contributions of the order g'^ and higher in and 7^ are not forbidden. 



5. Fixed points 



It is well known that IR asymptotic behaviour of a multiplicatively renormalizable field 
theory is governed by IR attractive fixed points of the corresponding RG equations. 
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Their coordinates are found from the requirement that all the /3 functions vanish, 
f^iid*) = 0) while the type of the point is determined by the matrix VL = {flik = 
d/3i/dgk\g=g^}: for an IR attractive fixed points it is positive, that is, the real parts of 
all its eigenvalues are positive. Here g = {gi} is the full set of couplings and /3j is the 
full set of the corresponding /3 functions. 

From the explicit expressions ( 14. 5 p and (14. 9 p for /3g it follows that the model (13. 2 p 
has a nontrivial fixed point 

9*S. = y^^)+0{f), (5.1) 

which is positive and IR attractive {dgf3g > 0) for ?/ > (of course, this fact is well 
known, see e.g. [2T], [22]). Substituting (15. ip into the equation (3u = and using the 
explicit expressions (14. 5 p and (14.91) gives 

2Q = u{u + iyd{d-l), (5.2) 

with Q from (14. 7p and corrections of order 0{y). 

The last equation is = 0. From eqns. (14. 5 p and (14.91) one finds that it is satisfied 
automatically up to the order 0{g). Thus there are two possibilities that cannot be 
distinguished within the one-loop approximation: 

The first one is that = 1 to all orders in g, as it happens in the rapid-change 
version of our model [Hj. Then the equation = becomes an identity and imposes 
no restriction on the coordinates of the fixed points. Then eq. (15.20 determines the 
coordinate function of the remaining free parameter A. 

In particular, for the most interesting physical case d = 3, the simple numerical 
analysis shows that the positive solution -u* of eq. ( 15. 2p is unique and exists for all 
A. As a function of A, it achieves a minimum u^, ~ 0.94 for A — —0.5 and grows 
as = 1^1 + 0(1) for A — i- ±oo. Some special values are u^: ~ 1.393 for ^ = 1 in 
agreement with the kinematic fixed point of the full-scale MHD problem [21], m* = 1 
for ^ = in agreement with [13] and m^, = 1 for A = —1. The simple inspection shows 
that this fixed point is IR attractive: duf3u > 0, 9^/3^ = = 0. 

A very similar behaviour of the solution u^,{A) takes place for all d > 2. As a 
function of d, it decreases monotonically and tends to unity as d tends to infinity. 
The explicit analytic solution of the cubic equation (15. 2p for general d looks rather 
cumbersome and we do not present it here. For d < 2 our results become inapplicable 
because the renormalization of the NS model ( 13.21) itself must be revisited [22] . 

Another possibility is that the function /3_4 has a nonvanishing contribution of 
order g'^ or higher. Then the equations Pu = Pa = ^ determine the possible values 
of the coordinates and A^. To find all their values, the two- loop calculation of Z_4 
is needed. However, it is clear without any calculation that = and A^ = 1 
are among the possible fixed-point values of A to all orders in g: the first case 
possesses additional symmetry with respect to the shift 9^9 + const (only derivatives 
of 9 enter the stochastic equation ( 12. ip ). while for the second case the nonlinearity 
Vi = {vkdk)9i — {9kdk)vi = dk{vk9i — 9kVi) in (12. ip is transverse: diVi = 0, so that 
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the nonlocal pressure term (12.31) vanishes. The both properties are preserved by the 
renormalization procedure. 

Existence of the fixed points different from A* = and 1 and their stability cannot 
be established without the explicit two-loop calculation. This issue lies beyound the 
scope of the present paper; here we only can say that for the passive vector field advected 
by the compressible Kraichnan's ensemble such points do exist; see [15]. 

6. Inertial-range anomalous scaling of the correlation functions, composite 
fields and operator product expansion 

The key role in the following is played by the critical dimensions A„ ^ associated with 
the irreducible tensor composite fields ( "local composite operators" in the field theoretic 
terminology) built solely of the fields ^ at a single space-time point x = (t, x). They 
have the forms 

F^^i = e,,{x)---e,Xx) {e,{x)e,{x)f + . . . , (e.i) 

where I < n is the number of the free vector indices and n = / -|- 2p is the total number 
of the fields 6 entering into the operator; the tensor indices and the argument x of the 
symbol Fn^i are omitted. The ellipsis stands for the appropriate subtractions involving 
the Kronecker delta symbols, which ensure that the resulting expressions are traceless 
with respect to contraction of any given pair of indices, for example, 9i9j — 5ij{6k6k/d) 
and so on. 

The quantities of interest are, in particular, the equal-time pair correlation functions 
of the operators ( 16. ip . For these, solving the corresponding RG equations gives the 
following asymptotic expression 

x)Ffcj(t, x')) - r-^-'-^'^- Cn,i;kMr) (6.2) 

with r = |x — x'l and certain scaling functions ({mr). To simplify the notation, here 
and below in similar expressions we omit the tensor indices and the labels of the scaling 
functions; the IR irrelevant parameters (like A or uq) are also not shown. 

The last expression in (16. 2p is valid for Ar ^ 1 and arbitrary values of mr. The 
inertial-convective range corresponds to the additional condition that mr ^ 1. The 
forms of the functions ({mr) are not determined by the RG equations themselves; their 
behavior for mr — )■ is studied using Wilson's OPE. 

According to the OPE, the equal-time product Fi{x)F2{x') of two renormalized 
composite operators at x = (x -|- x')/2 = const and r = x — x' — >■ can be represented 
in the form 

F,{x)F2{x') ~ (r)F(t,x), (6.3) 

F 

where the functions Cp are the Wilson coefficients, regular in m^, and F are, in 
general, all possible renormalized local composite operators allowed by symmetry. 
More precisely, the operators entering the OPE are those which appear in the 
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corresponding Taylor expansions, and also all possible operators that admix to them in 
renormalization. If these operators have additional vector indices, they are contracted 
with the corresponding indices of the coefficients Cp. 

It can always be assumed that the expansion in Eq. (16. 3 p is made in operators 
with definite critical dimensions Ap. The correlation functions fl6.2p are obtained by 
averaging equation of the type (16.30 with the weight exp5, where S is the action 
functional (14. 2p : the quantities (F) appear on the right hand sides. Their asymptotic 
behavior for m — )■ is found from the corresponding RG equations and has the form 
(F) oc m^-p. 

From the expansion (16. 3 p we therefore find the following asymptotic expression for 
the scaling function ({mr) in the representation (16. 2p for mr ^ 1: 

C(mr) ~ ^ (mr)^^, (6.4) 

F 

where the coefficients Ap = Ap{mr) are regular in {mr)'^. 



7. Anomalous scaling and the exponents in the one-loop approximation 

The feature specific of the models of turbulence is the existence of composite operators 
with negative critical dimensions. Such operators are termed "dangerous," because their 
contributions to the OPE diverge at mr — )■ [211 122] . 

Obviously, most dangerous are the operators (16. ip with the critical dimensions 
A.n,i = —n + 0{y). Like in the Kraichnan's case, the analysis shows that their anomalous 
dimensions can be calculated in the simplified model without the random forcing in 
the stochastic equation (12. ip because the correlator (12. 2p does not enter the relevant 
Feynman diagrams [9]. Then those operators become multiplicatively renormalizable. 



Fn^i = Zn,iF^i. The practical one- loop calculation of the renormalization constants Zn,i 
is similar to the case of Kraichnan's velocity field, discussed in [9l [15] in detail, so that 
here we give only the result: 

gSd A'Qni 1 



u{u+ 1) 4:d{d + 2) y 



0(9'), (7.1) 



where 



Qn,i = 2n{n - 1) - (d + l){n - l){d + n + I - 2) = 

= -{d-l){n-l){d + n + l) + 2l{l - 1) (7.2) 

and Sd is defined below equation (14. 7p . Note that the same polynomial Qn,i arises in 
the scalar case |25] and in Kraichnan's MHD model [9]. 
The corresponding anomalous dimension is 



^-'^ = -f^\ i:77;7Tk + ^^3'). (7.3) 

u[u + 1) Ad[d + 2) 
Substituting the fixed-point value ( 15. ip gives 

u{u + 1) Sd[d — 1) 
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If the function vanishes identically, the solution u^{A) of the equation (15.21] 
should be substituted into (17.41) : then the dependence on the free parameter A persists 
in 7„^/. Otherwise, the fixed-point values m*, A^, should be used; see discussion in 
section [5l In particular, for A* = I the fixed-point value of is the positive solution 
of the quadratic equation u{u + 1) = 2{d + 2) / d. Then eq. (17.41) becomes 

^"''= 6(ci-?)7^ + 2) ^ + ^^^'^' ^^-'^ 

which agrees with the result derived in [19j for the MHD case. Note that (17. 5p coincides 
with its analog in the Kraichnan's case [9] up to the substitution ^ — )■ y/3. 

For ^ = 0, the anomalous dimensions 7„^i vanish to all orders in y, because the 
operators F„ ; become UV finite and are not renormalized. In that case, interesting 
quantities are structure functions (rather than plain correlation functions); their inertial- 
range behaviour is determined by the operators built of the derivatives of the fields 6; 
see [121113]. 

With this exception, the amplitude A^/u{u + 1) in (17. 4p is positive for any physical 
fixed point. Thus the dimension is negative for the most interesting case of the 
scalar operator with I = and increases monotonically with / (for a fixed n). 

From the relation A„ ^ = —n + 0{y) it follows that the critical dimensions satisfy 
the same hierarchy relations: An,i > A„_// if / > which are conveniently expressed as 
inequalities dAn,i/dl > 0. 

This fact, first established in |8] for the Kraichnan's MHD model, has a deep 
physical meaning: in the presence of large-scale anisotropy, the leading contribution 
in the inertial-range behavior mr — > of the correlation function like (16. 2 p is given by 
the isotropic "shell" (/ = 0). The corresponding anomalous exponent is the same as 
for the purely isotropic case. The anisotropic contributions give only corrections which 
vanish for mr — )■ 0, the faster the higher the degree of anisotropy / is. This effect gives 
some quantitative support for Kolmogorov's hypothesis of the local isotropy restoration 
and appears rather robust, being observed for the real fluid turbulence [26] and the 
passive scalar model [25] . 

8. Conclusion 

We have studied a model of a divergence-free (transverse) vector quantity 6, passively 
advected by a random non-Gaussian velocity field with finite correlation time, governed 
by the stochastic NS equation. The model is described by an advection-diffusion 
equation with a random large-scale stirring force, nonlocal pressure term and the most 
general form of the inertial nonlinearity, "controlled" by the parameter A oc ^o- 

We have shown that, in the inertial range of scales, the correlation functions of the 
field 6 exhibit anomalous scaling behaviour. The corresponding anomalous exponents 
are determined by the critical dimensions of tensor composite fields (16. ip built solely of 
the passive vector field. They are calculated (including the anisotropic sectors) to the 
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leading order of the expansion in y, the exponent entering the correlation function of the 
stirring force in the NS equation (the one- loop approximation in the RG terminology). 

Like in the special MHD case A = Aq = 1, the exponents exhibit a kind of hierarchy 
related to the degree of anisotropy: the less is the rank of the tensor operator, the less is 
the dimension and, consequently, the more important is the contribution to the inertial- 
range behaviour. Thus in the presence of large-scale anisotropy the leading terms, 
determined by the scalar operators, remain the same as in the purely isotropic case, in 
agreement with the phenomenological hypothesis of the local isotropy restoration. 

The question that remains open is whether the amplitude A in front of the 
"stretching term" {9d)v in the advection-diffusion equation tends to some fixed-point 
values, or it remains a free parameter which the anomalous exponents depend upon. 
The analysis of that alternative lies beyound the scope of the one-loop approximation. 
We plan to perform it in the nearest future. 
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